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Señales y sistemas : Apuntes procesado de señal 


1. Señales y sistemas. 


En cualquier transmisión de información hay dos elementos fundamentales: las 
señales, que serán las que contendrán la información y los sistemas, que son los encargados 
de manipular las señales y los que las conducirán de la fuente al destino. En este primer 
tema se realizarán toda una serie de definiciones, tanto sobre señales, como sobre sistemas, 
que permitirá caracterizarlos. 


1.1. Señales. Tipos. 


Hoy en día, cuando se oye hablar de señales uno automáticamente piensa en señales 
de naturaleza electromagnética, durante estos apuntes, salvo contadas excepciones, no se 
asignarán unidades a las señales, se las tratará de forma puramente matemática. Se realiza 
esta advertencia porque alguna de las definiciones que se realizarán dentro de este apartado, 
puede que induzca a pensar en un significado físico, sin embargo, no debe darse nunca un 
sentido físico. 

Respecto a la naturaleza de las señales, se pueden clasificar de dos formas 
diferentes: 

Una primera clasificación será entre señales deterministas y aleatorias: 


Deterministas: son aquellas señales que admiten una expresión matemática que las define 
para cualquier instante de tiempo. Este tipo de señales no aparecen nunca en una 
comunicación, no es lógico transmitir una señal de la que se conoce su expresión 
matemática, sería más rápido trasmitir la función matemática de la define. Las señales 
determinista serán utilizadas para testear los sistemas y conocer si estos funcionan como es 
debido, o por el contrario, se utilizarán para caracterizar sistemas no conocidos. 


Aleatorias: son señales que no atiende a ninguna expresión matemática, sino que en cada 
instante de tiempo toman un valor, entre un conjunto de valores determinado, con una cierta 
probabilidad. Este tipo de señales constituyen las comunicaciones reales. Los parámetros 
que se puede conocer de las señales aleatorias son estadísticos: media, varianza, densidad de 
probabilidad, etc. Como se puede suponer, estos datos no son suficientes para poder 
comprobar si un sistema funciona o para caracterizar un sistema. De ahí que sean 
importantes las señales deterministas. j 

La segunda clasificación será entre señales analógicas (continuas) y discretas, 
digitales o numéricas. Dicha clasificación es doble, se puede realizar respecto a la variable 
tiempo o respecto al valor que va tomando la señal en los distintos instantes de tiempo. Las 
señales cuya discretización es sobre la variable tiempo suelen llamarse secuencias. 


Analógica: son aquellas señales que entre un valor máximo y un valor mínimo tienen 
infinitos valores intermedios. 


Discreta: son aquellas a las que entre un valor máximo y uno mínimo sólo les está 
permitido tomar un número finito de valores. 
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Señal continua en tiempo y discreta en 
valor. 


Señal continua en tiempo y en valor. 


x(t) 


x(0) 


continua en valor, 
x(n) aj 


1.1.1. Valores característicos de las señales. 


Las señales, independientemente de como estén enmarcadas dentro de las 
clasificaciones que se acaban de realizar, se las caracteriza por diversos parámetros. Los 
más usuales son: 


Valor instantáneo: es el valor que toma la señal en cada instante de tiempo. 


x(£) 


Valor de pico: es el valor máximo, tanto positivo como negativo, de la señal. Se suele 
expresar siempre como un número positivo. 
O 


Valor medio de la señal: es el valor que tomaría la señal si sus valores se repartieran por 
igual en todo el eje de tiempos. Tiene sentido en señales infinitas en tiempo. 


1 fra 
Jim 7 le x(t) dt 
En señales periódicas, la expresión anterior queda reducida a un periodo: yl 
7 pe x(1) dt donde 7 es el periodo de la señal 
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Potencia instantánea: a este parámetro no hay que buscarle un significado fisico, no lo 
tiene, pese a que su nombre puede indicar lo contrario. Está definida como: 


Oj 


Potencia de pico: es la potencia instantánea máxima. 
2 
OI 


Potencia media: tiene el mismo significado que el valor medio de la señal, pero con la 
potencia. Su definición es: 


1 
pro a 


En señales periódicas, la expresión anterior queda reducida a un periodo: 
1 fr/ 


7 O di donde T es el periodo de la señ al 


Energía: tiene sentido en señales finitas en tiempo. Al igual que la potencia tampoco tiene 
significado físico. Su definición es: 


Plxco? as 


1.1.2. Señales más usuales. Definición. 


Señales analógicas en tiempo Señales discretas en tiempo 
La variable tiempo es una variable La variable n es una variable discreta, sólo 
continua, puede tomar cualquier valor. puede tomar valores enteros. 
Función escalón analógica: Secuencia escalón: 
u(t) u(n) 
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Función signo analógica: 


sgn(t) 


1 1>0 
t= 
sent) -1 t<0 
Función pulso analógica: 


TI(ut,) 


4/2 172 t 


1 ñ 
ns) 1 (SIS tp + 


h 0 resto 


Función triángulo analógica: 


A(us,) 


pliz) 1 SISt, +, Jn 
, Fl A ñ; 


Apuntes procesado de señal 


Secuencia signo: 


sgn(n) 


1 n>0 
sen(m)=)_, 


Secuencia pulso: 


IN) 


-n, /2 n /2 n 
UN 


A 
nz). 1 M7 SNS + 2 


A O resto 


Secuencia triángulo: 


hr» 
a NN SNSN/+I 


0 resto 
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Función sinc analógica: - Secuencia sine: 
sinc(t) sinc(n) 
e. 0 
1 et 
t 6) ... e. en 
sinc(af) = er) sinc(an) = Lia 
art amn 
Delta de Dirac (impulso): Delta de Kronecker (Impulso discreto): 
90) 
Sn) 
A E E A 


Se puede definir de varias formas, entre 
ellas se pueden destacar: 1 n=0 
s(n)= 


lim8,(0) = Xt) 0 resto 


La Delta es un caso especial, es lo que en términos de álgebra se denomina un 
funcional, es decir, una función que sólo aparece bajo el signo de una operación. En el caso 
de la Delta de Dirac la operación bajo la que suele aparecer es la integral. Dada esta 
particularidad la Delta de Dirac se suele definir con estas características: 


Psmar=1  Fxost=todt=x(t5) 


En el caso discreto sucede lo mismo, con una pequeña diferencia; suele aparecer 
bajo el signo sumatorio y viene caracterizada por las propiedades: 


Y dm) =1 Y x(mAn=n,) = x(m) 


n=-00 n=-0 


1.2. Sistemas. Tipos. 


Se puede definir un sistema como un conjunto de elementos que tienen una 
determinada finalidad, pero desde el punto de vista que interesa en procesado de señal, 
también se puede definir como una ley de correspondencia que asigna a un conjunto de 
entradas, x;(t), un conjunto de salidas, y;(t). Por tanto, se puede caracterizar un sistema por 
una función matemática, S[+], tal que: 
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y() =S[x(0)] 


Los sistemas, dependiendo de las señales que procesan, se pueden dividir en 
sistemas discretos o digitales y sistemas analógicos. 


1.2.1. Características de los sistemas 


Linealidad. : e 
S|ox,(£) + Px, (9] SS aS|x, (2)] + 8S[x, (1)] 


Esta característica es muy importante en procesado de señal, pues, permitirá 
encontrar la salida de un sistema a una entrada complicada, como la suma de salidas a 
entradas más sencillas. 


Invarianza temporal. 
Slx(1-7)]=y(-7)  paraWVreRe 


Los sistemas que tengan esta característica más la de linealidad podrán ser 
caracterizados a partir de una señal de entrada y una señal de salida. 


Causalidad. (temporal) 
Six(1)=0 parat<t, >y(1)=0 Vi<t, 


En realidad todos los sistemas físicos son causales, pues no es posible responder a 
una entrada antes de tener entrada. Sin embargo, matemáticamente es posible definirlos, por 
tanto, es algo que se tendrá que tener en cuenta al diseñarlos. 


Estabilidad. (dentro de sistemas lineales) 
Silx(1))<M<o=>|y(0)]<M' 


Invertibilidad. 
Un sistema será invertible sí, definiendo el sistema como y(t)=S[x(t)], existe otro 
sistema tal que S;¡[y(t)]=k*x(t-tp). 


ás LA Etuta 


Existen más características propias de los sistemas, pero las nombradas 
anteriormente son las más importantes. Las características enunciadas se definen por igual 
en sistemas discretos y en sistemas analógicos. 
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1.2.2, Respuesta impulsional. 


Se define como respuesta impulsional de un sistema la respuesta del sistema cuando 
la entrada del sistema es una Delta de Dirac. 
Si x(t)=0(t): 


h(t)=S[S(t)] Respuesta impulsional 


Si el sistema no es invariante en el tiempo: h(t,tp)=S[S(t-tp)] 
El sistema, en este caso, tendrá una respuesta impulsional diferente para cada 
retardo tg. 


Aunque la respuesta impulsional, por definición, se podría encontrar para cualquier 
sistema, en el siguiente apartado se mostrará que sólo tiene sentido definirla en sistemas 
lineales, por tanto, siempre que se hable de. respuesta impulsional se entenderá como 
perteneciente a un sistema lineal. 


1.3. ' Sistemas lineales e invariantes en el tiempo. (L.T.I) 


Los sistemas L.T.I. permiten ser caracterizados por completo por la respuesta 
impulsional. 
Si se dispone de un sistema L.T.I., al que se excita con una señal x(t): 


0 


yo =stx0]=s[f 50 -0d] = [ sixcyst- ola: = PP xoslsa le = 
E 1 
x(1)= IF x(m)S(1-1)dt linealidad 


= pe x(t)h(t,tjdt 


Como además, el sistema es invariante en el tiempo: 
h(t,7) =h(t- 7) 


por tanto, se obtiene como resultado final: 
y) = Pon - wat 


Esta expresión es conocida como integral de convolución. 

Haciendo un resumen del apartado, en un sistema L.T.I., se puede conocer la 
respuesta a cualquier excitación, simplemente conociendo la respuesta impulsional, que 
caracteriza por completo este tipo de sistemas, y la excitación. Si el sistema no tiene las 
características L.T.L., la respuesta impulsional no caracteriza al sistema. 

En general, se suele utilizar la nomenclatura respuesta impulsional para la respuesta 
a una delta, solamente en sistemas L.T.L. 
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1.4. Convolución. 


La integral de convolución, normalmemte llamada convolución, está definida como: 
yo = E xo a 


Por simplificar la expresión la notación que generalmente se utiliza es: 
(0) =x(0*h(0) 


1.4.1. Propiedades de la convolución. 


Conmutativa. 


y (0) =(0*h(1) = h(0)* x(1) 
y=Fxona-ya =$ noxe 


Para demostrar esta propiedad simplemente hay que realizar el cambio de 
variable: 7=1-— 7 


Asociativa. 


YO) = Ay ya =x0* 000) 


Para demostrar la propiedad asociativa, se deben sustituir las convoluciones por su 
expresión integral y cambiar de orden las integrales. 


Distributiva. 


(aux(2) + By(n)ez(1) =ax(1)z(1) + By(o)*z(0) 

Todas las propiedades que aparecen en esta relación y en particular la distributiva, 
son fácilmente demostrables si se tiene en cuenta que la convolución no es más que una 
integral y, por tanto, tendrá las propiedades de la operación integral. 

Derivación. 


ON O) Lo) YU) = xo Lo) 


Elemento neutro. 


x(0)* (0) = 1) x(0) = x(0) 


Esta propiedad era una de las características de la Delta de Dirac. 
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1.4.2. Ejemplo: Resolución gráfica de la convolución. 


Calcular la convolución de: 


x(0) h(t) 
: t t 


(2 


x(1)= HA h(t)= rn 


yo= E xoni-oa 


Primer paso: 

e Realizar el cambio de variable T=—T para invertir una de las señales, normalmente 
la que sea más fácil, recuerde que la convolución tiene la propiedad conmutativa. 
Gráficamente representa la inversión de sentido de una de las funciones. 


x(r) h(-0) 
T T : 
Segundo paso: 


Desplazar la señal que se ha invertido desde -co hasta oo, para, posteriormente, ir 
calculando la integral. Para forzar el desplazamiento de forma fácil, es conveniente realizar 
el cambio de variable: -7=1- r, en la función que en el primer paso se ha invertido. 


Tercer paso: 

| El último paso es calcular la integral del producto de las dos funciones para cada 
| desplazamiento t. . 

| A la hora de realizar el cálculo, conviene realizarlo por tramos, en los cuales los 
| límites de la integral sean los mismos. 
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Aplicando este consejo al ejemplo propuesto, el tercer paso, se subdividirá en tres 
tramos. 


Primer tramo. 


t<0; no hay solapamiento de las dos funciones. 
La integral valdrá siempre cero, 
y(t)=0- para t<0. 


Segundo tramo. 


O <t< T; hay un solapamiento parcial, uno de los límites de la integral es siempre 
cero y el otro depende del e t. 


il 
== “dí= 2 1-e/% 0<t< 
MO) liz dt pu e T 


Tercer tramo: 


> T; hay un solapamiento en todo el pulso, los dos limites de la integral son 


constantes. 
1] --» Ét T, = 
»(M= lie dr = “e - 1)e %  t>T 
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1.43. Ejemplo propuesto. 


Calcular: 
(1) =x(0*y(0) 
Í 1 
x(1)= 414] y) = 411 — T, >1) 
Y T 
Resultado: 

z(t) 
Cc 

A -B - B A t 

Ti E ET 
sd HA E eE ; C=AA,7, 


1.5. Propiedades de los sistemas a partir de la respuesta impulsional. 


Una vez definida la respuesta impulsional, las propiedades de los sistemas y de la 
convolución, hay tres propiedades de los sistemas que pueden reconocerse a partir de la 
respuesta impulsional. Queda presupuesto que en un sistema, para que exista respuesta 
impulsional, entendida como tal, el sistema tiene que ser lineal, 


Invarianza temporal. 
Un sistema será invariante en el tiempo si su respuesta impulsional sólo depende de 


la variable tiempo, no del instante de tiempo para el cual se defina, es decir, h(t,tp)=h(t-tp). 


Causalidad. : 
Si se aplica la definición de causalidad al caso particular de una excitación S(t), el 
sistema será causal si, h(t)=0 parat<0, 


Estabilidad. 
Si ]x(1)] <M 


b|= E sC0hti> var SÍ Iitolecr 0) lee < MÍ” coles 


Para que y(t) este acotada: 
Plhoje <o0 


La condición encontrada es suficiente y necesaria, 
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2. Análisis de Fourier. 


En el capítulo anterior se ha visto que si un sistema era lineal se podía calcular su 
respuesta a una excitación con una expresión matemática complicada, como suma de 
respuestas a señales que fueran más fáciles de tratar, En este capítulo se estudiará una forma 
sistemática de realizar esta descomposición, de forma que sea fácil tratar las señales y por 
supuesto los sistemas. 


2.1. Series de Fourier. 


Los Desarrollos en Serie de Fourier, no son más que la descomposición de una 
señal, que cumpla una serie de condiciones, en un sumatorio de señales ortogonales. Se 
denominan señales ortogonales aquellas que cumplen: 


«(1 y" (Ddt=0 excepto si y(t) = Kx(1) 


En el caso de los Desarrollos en Serie de Fourier, la familia de señales ortogonales 
que se utilizarán serán: exponenciales complejas, o sí se prefiere, senos y cosenos. 

En general, para hacer una descomposición, no es necesario que las funciones 
ortogonales sean siempre exponenciales complejas, puede ser cualquier familia de funciones 
ortogonales. Más adelante se comprenderá la razón por la que se utilizan este tipo de 
funciones y no cualquier otra. 

Las condiciones que tiene que cumplir una función matemática para poder ser 
aproximada por una Serie de Fourier, son: 


1) Ser periódica de periodo 7, = ello E 
fo % 


x(1 +mT,) = x(£), para V t, siendo m entero 


2) La función tiene que estar definida para todo t. 
=0<tf<0w 


3) Ser continua en [o, T,] o tener un número finito de discontinuidades de primera especie. 


EXP ES | 


E e E 
Discontinudad de primera especie Discontinuidad de segunda especi 
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4) Tener un número finito de máximos y mínimos en [0, T,]. 


Si se cumplen todas estas condiciones se puede demostrar-(por aproximación de 
funciones por mínimos cuadrados, con funciones ortogonales) que: 


x(£) = + da, cos(n0,1)+b, sen(70,£)) 
n=1 


o utilizando las fórmulas de Euler: 
Formulas de Euler => eo - cos(wyt) + jsen( opt) 


de forma más compacta: 


E : 
x(1) = S ceros 0 =27% 


n=-0 


Esto es lo se denomina Desarrollo en Serie de F ourier, 
En procesado de señal, cada uno de los términos del sumatorio, recibe un nombre 
especial, dependiendo del valor de n. . 
El termino n=0, se llama valor medio o continua. 
El termino n=1, se llama fundamental. 
Los demás términos se llaman armónicos. 
Los coeficientes del desarrollo se encontrarán con la expresión: 


1 -jno, 
C,= ql x0e 0 e 


En general, los coeficientes Cy, serán complejos, C,, =|C, e? ae[C,] 

El cálculo de los coeficientes del Desarrollo en Serie de Fourier admite bastantes 
simplificación, si se tiene en cuenta las características típicas de las funciones matemáticas 
que se utilizan en procesado para representar las señales. 

-Si x(t) es real: 4 


1 si 
Eo 0er C,= [ xoe “dt 


“Y 
7, 


como x(t) es real: 


als 


arg[C,] = -arg[C.,] 


6,=0,>| 


Si se cumplen estas igualdades, se dice que tiene simetría Hermítica. 
Los coeficientes se suelen representar mediante dos gráficas, una para la fase 
(espectro de fase) y otra para el módulo (espectro de amplitud). 


AA A 
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A 4 0 £ 24 ND arg[C,] 


En las gráficas se puede observar que la representación, tanto del espectro de fase 
como del de amplitud, sólo tiene valores distintos de cero en determinados valores de 
frecuencia, cada uno de esos valores se le llama raya espectral. También se puede apreciar la 
situación de estas rayas espectrales, están equiespaciadas y situadas a múltiplos de la 
frecuencia fundamental fo: O, £o, 20, 3£o, ..., lo que se ha denominado como armónicos. 

Si x(t) es par (x(t)=x(-t)): 

Los coeficiente Cy serán reales y, por tanto, su fase será 0 ó 180". 
La expresión que permite calcularlos se simplifica, siendo: 


1 
C, = Te x(1) * cos(n0,1)dt 


Si x(t) es impar (x(t)=-x(-t)): 
Los coeficiente Cy serán imaginarios puros y, por tanto, su fase será 90% ó -90%, 
La expresión que permite calcularlos se simplifica, siendo: 


e . 
C, = T, [ xo) sen(no,t)dt 


Para demostrar estas dos simplificaciones simplemente hay que utilizar las fórmulas 
de Euler, y aprovechar las simetrías del seno (impar) y el coseno (par). 


2.1.1. Ejemplo: Cálculo del Desarrollo en Serie de Fourier . 


x(t) 


-T, /2 -112 2 1,2 t 


2 de =h _1 1 JN e -jnoy/2 _ ¿jnoyt/2) — 
Cj= T, [ xo Jo dy = 7 Ese Fo de = le e ) 


A 
- Tono, 
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Utilizando las fórmulas de Euler: 
eh 


2) 


= A sen(1mf,7) 
mm 


sen(x) = 


Resultados de este tipo se pueden compactar si se tiene en cuenta la definición de la 
función sinc: 
sen( e 
sinc(x)= Lio 


7 
C, = Afporsinc(nfyT) 


Si se representa gráficamente se tiene: 


A mod[G,] 
py! 


Por tanto, x(t) se puede escribir como: 


Ñ x(t)= eN Aforsinc(mf) 7) eme 


n=-00 


A A 
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2.1.2, Ejemplo propuesto. 


Encontrar el Desarrollo en Serie de Fourier de: 


x(t) 
A 
X, /4 7/4 [37 /4 
Ñ t 
-A 
a+b a-b 2 
Nota: e“ te? =e 2107 40” al : 
Resultado: 
fo ,»n=0 
Lac 
El= 4sinb?) ,N impar 
lo ,n par 


x(1) = Cy + D /2C, cos(27mfot) 


n=1 


2.2. Teorema de Parseval. 
Por definición, la potencia media de una señal periódica es: 


1 2 
==, lx(oP ar 


Desarrollando esta igualdad, teniendo en cuenta que la señal sobre la que se realiza 
el cálculo es periódica: 


P= AL eor dt= Lxox 0d ==, 0) Ecprma| dt= 
=2f,x0] Eco roo la - $ 27 - x(0)C; "2 dp 


e do 


= Eq + x(0e 2 a - E 
n= o 


n=-0 


En definitiva, el teorema de Parseval proporciona una segunda forma de calcular la 
potencia media de una señal: . 
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P=—Í hwPa= lc] 


L, no 
Ejemplo: Cálculo de la potencia media de una señal. 


C, 


2.2.1, 


Dados los espectros de fase y de amplitud de las figuras, calcular la expresión de la 


señal temporal correspondiente (periódica). Calcular la potencia media en los dos dominios, 
temporal y frecuencial. 


100 200 f(Hz) 


-200  -100 0 100 200 RHz) . 2 


00 
x(1)= Dyer y =2x% 


n=-00 
f)=100Hz 

Co =-] 

C al C; =2 

C=j 

Co=-j 


x(1) =-1420/2%8 420700 y ¡e 24 - jose 5 
=-142(0/200 4) y ¡(OSA ¿RA 
x(t) = -14+2c0s(2 11001) —-2 sen(272001) 


Cálculo de la potencia media en el tiempo: 


1 2 
=xÍ, kof ar 


Al ser un sumatorio de armónicos, la potencia media es la suma de potencias de cada 
uno de los armónicos (ortogonalidad). 


p* cos(mo,t) sen(no,t) dt =0 V m,n 
f a cos(mo,t) cos(10,t) dt =0 Vmin 
f sen(mo,f)sen(no,1)dt=0 . Vmzn 
» 7 T, 
f sen*(no,£) dt = f cos? (10,£) dt = E 
a 
P(-D)=1 lim Ly dt=1 


A 5 
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1 1.7% 
P(4cos(2111001)) =8 16)” cos (211000 di = 162% 
Ls T 2 


Ta 
4 


lo fm 
P(2sen(212001)) = 2 ql, sen* (272001) dt = E 
0 


ES 
1, 
P=1+8+2=1lw 


Cálculo de la potencia media sobre el espectro: 


00 
P= Y [C,' =1+22+2?+141=11w 


n=-0 


Evidentemente, el resultado es idéntico para los dos casos. 

El teorema de Parseval no siempre es la mejor forma de calcular la potencia media 
de la señal, sino que simplemente es otra opción, depende de la señal se debe utilizar un 
método u otro, obiamente la que presente menor dificultad matemática. 


2.3. Respuesta en régimen permanente sinusoidal. 


Se ha visto como descomponer una señal periódica en un sumatorio de 
exponenciales complejas. En este apartado va a darse una justificación a este tipo de 
desarrollos. 

Primeramente se verá como reacciona un sistema L.T.I. a una excitación del tipo 


exponencial compleja. 


x(1) = e/% 


La salida y(t) se encontrará realizando la convolución de la respuesta impulsional del 
sistema con la entrada: 


yo = 00 =Í xaomoar= [ema en lnea: 
=e"H(jo,) HGo,)= Prgera =cte 


De este cálculo se pueden deducir una serie de particularidades de los sistemas L.T.I 
y las exponenciales complejas, también llamadas fasores. 

- La exponencial compleja es una autofunción de la convolución. Una autofunción es 
una función que sobrevive a una operación, en este caso se tiene que después de realizar la 
convolución se tiene la misma exponencial que a la entrada, solamente afectada por un 
cambio de amplitud (Hop). 

- Por otra parte, un sistema L.T.I. reacciona ante un fasor, simplemente cambiando 
la amplitud de este (amplitud que puede ser compleja, módulo y fase). En general, si se 
dispone de una entrada de la forma x(f) = e“”, donde w puede ser cualquier pulsación, y se 
realiza la operación: 

IÓ rin 
0 H(jo) 
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la expresión H(¡o) representa como responde el sistema a cada pulsación en particular, a 
esta expresión se le llama función de transferencia del sistema. Por simplificar la notación, la 
función de transferencia del sistema, también llamada respuesta en frecuencia del sistema, se 
suele expresar como, H(f). 

-Si se dispone de un sistema L.T.I, que se excita con una señal periódica, la cual 
admita un desarrollo en serie de la siguiente forma: 


x(9= Y Cero a =2%, 


n=-0 


y se quiere conocer la salida, y(t), como el sistema es lineal se puede calcular la salida como 
suma de las respuestas del sistema a los diferentes fasores que componen la entrada: 


10-38 Eso ¡exo= Elstxco) 


en este caso particular: . 


10=3 Ec. [o ym = Zls[c,e””"!) 


n=-00 


Como el sistema es lineal e invariante, la respuesta a cada una de las excitaciones se 
puede calcular como la convolución de la respuesta impulsional del sistema y el fasor 
correspondiente. La repuesta a un fasor ha sido calculada al principio de este apartado y 
era: 


S[e]=e"HGjo,)  Hljo,)= Panor =cle 


Luego en el caso de una excitación periódica que admita desarrollo en serie de 
Fourier, la respuesta de un sistema L.T.I. se podrá calcular como: 


yo=8 Ec Jero= Ec mue 


n=-0 


Por tanto, la respuesta de una señal periódica, o en general la respuesta a una señal 
compuesta de un sumatorio de fasores (el sumatorio de una serie de señales periódicas es 
una señal periódica), puede ser calculada como el sumatorio de los mismos fasores, sólo 
que cada uno de ellos afectado por una constante, H(fp), propia del sistema. Si el sistema no 
es L.T.I. este desarrollo no es cierto. 

- Si el desarrollo anterior se aplica al cálculo de la potencia media de una señal, en 
general, compuesta de un sumatorio de fasores, la respuesta del sistema, será calculada 
como: +. 


y) = DY CH) 


n=-0 


El teorema de Parseval decía que si se disponía de una señal periódica, se podía 
calcular la potencia media como: 


e Ec) > 


$ 2 


no 


C, 


n 
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En este caso los coeficientes Cyn de la señal y(t), serían C,H (nf), por tanto, la 
potencia media de la señal y(t) sería: 


Br) = Y lc, Haro = DIC ar 


n=-0 n=-00 


De este apartado se puede destacar la operación: 
AGO) = mg ra 


que se conoce como Transformada de Fourier. 


2.4. Transformada de Fourier. 


La Transformada de Fourier, X(f), de una señal de energía finita, extensible a señales 
de potencia media finita, (si se admite la existencia de Deltas de Dirac), está definida como: 


F(x(0) = X()= Puerta 


Puede verse que no es más que una Transformada de Laplace, bilateral, 
particularizada para S=jw. El hecho de definirse como bilateral, no es más que permitir que 
la señal esté definida desde -co y hasta oo, es decir, no se tendrán en cuenta ninguna 
condición inicial, al contrario que sucedía en la Transformada de Laplace. 

A su vez puede calcularse la Transformada Inversa de Fourier: 


FX) =x0=Í Xd, donde X(7)=FÍx(0)) 


Tanto la transformada, como la transformada inversa son relaciones biunívocas, es 
decir, que no hay dos señales que tengan la misma transformada, ni dos señales en dominio 
transformado que tengan la misma transformada inversa. 

Una justificación del cálculo de esta transformada se puede buscar en las Series de 
Fourier. Una señal periódica, permitía un Desarrollo en Serie de Fourier de la siguiente 
forma: , 


xX0= DC, ; 0=2% ; C,= —Í, (1er dy 
s pa=eo' y» 
Se vio que las señales periódicas tenían rayas espectrales cada fp (fp era el inverso 
del periodo de la señal) y que cada raya espectral tenía amplitud C¡. Evidentemente, los 
Desarrollos en Serie de Fourier sólo son válidos para señales periódicas, pero si se dispone 
de una señal no periódica, se puede considerar que la señal tiene periodo infinito. Si se 
aplica esta consideración a un desarrollo en serie, como Ty es infinito, fp es cero, por tanto, 
las rayas espectrales están infinitamente juntas. Dicho de otra manera, no se tiene una serie 
de coeficientes Cy, sino una función continua C(f). Si se aplica esto las expresiones de los 
Desarrollos en Serie de Fourier, se tendrá: j 


21 , 
x(t)= 27 l ; x(1) + eo gr. giros 
alo 


pa 
C, . 
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Si ¿ >0<> f >0, los coeficientes del Desarrollo en Serie de Fourier pasan a ser una 
función continua: 


12 , z 
x(1)= lim s/ (me "grero" 
T,>w0 y pa) 
cl) 


Queda un sumatorio de términos infinitesimales, lo cual, no es más que la definición 
de integral: 


x(t) = E pl «(e dejes ay 
CUFXU) 


X(£) juega el papel de los coeficientes Cy, de los Desarrollos en Serie de Fourier. 

En señales periódicas: espectro discontinuo —> sumatorio. 

En señales no periódicas: espectro continuo > integral. 

Una demostración más rigurosa de las expresión de la transformada inversa se puede 
realizar basandose en la siguiente igualdad: 


1 o Si 
pl Jot 
51) 7 F. do 


En efecto: 
l fo, ol fa, . sen(af) ,  sen(at) 
a Jo = ap qe = A SO 
Lerdo Im>Fedo= m2 = lima O) 


“Para comprender mejor esta última igualdad se puede ver que la función, 
.. sen(at) 
lim a 

rat 


aso 


, tiene toda su área concentrada en el origen de coordenadas y vale la unidad: 


ota 


Si se tiene en cuenta esta igualdad: 
'o o. l fo, l fo, ffo , 
x«0)=Í x0sU-0ut= Ex Free Ja ss LE x(De "dido = 


1 00 y 
x(t)= La X(oJe""do 


De la misma forma que en las series, existen unas simetrías si x(t) cumple las 
siguientes condiciones: 
-Six(t) es real: 


Xxf =$ (De “d Xx n= nera 


como x(t) es real: 
IA 


X= Y E 
it oda =-argl Xf)] 


Una función que tenga estas características se dice que tiene simetría Hermítica. 
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La transformada se suele representar mediante dos gráficas, una para la fase 
(espectro de fase) y otra para el módulo (espectro de amplitud), donde se puede apreciar el 
significado de la simetria Hermítica. 
mod[X(0] 


— 


arg[X(D] 


— 


-Si x(t) es par (x(t)=x(-t)): 
La transformada será real y, por tanto, su fase será 0? ó 1800. 
La expresión que permite calcularlos se simplifica, siendo: 


Xx) =2L 0) cosa 


Si x(t) es impar (x(t)=-x(-t)): 
La transformada será imaginaria pura y, por tanto, su fase será 90% ó -902, 
La expresión que permite calcularlos se simplifica, siendo: 


X= O) sen(21/1)dt 


Para demostrar estas dos simplificaciones, simplemente hay que utilizar las fórmulas 
de Euler, y aprovechar las simetrías del seno (impar) y el coseno (par). 

Otra característica que presentan las señales en el dominio transformado es que si se 
particuraliza para £=0, el valor obtenido es el área de la señal de la cual se ha partido para 
obtener la transformada. Dicha característica permite solucionar integrales definidas cuyo 
integrando no tiene primitiva y sólo se podían resolver, aproximadamente, por métodos de 
cálculo numérico, por ejemplo, área de la función sinc. 


Far= Exnerar =X0 
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2.4.1. Ejemplo: Cálculo de la Transformada de Fourier de un pulso. 
X(0 


A 


4,12 t,/2 t> 


XK y )= f xd = 8 Aci" dt = ze =I% gim) psentaño) 
X(f)= Atysinc( ft) 


Como la señal es real y par, se puede calcular también de la siguiente forma: 


X(p= 170 cos(0t)dt = lu Acos(ot)dt = 24 ¿¿[2%). = At sinc(ft,) 


2.5. Teorema de Rayleigh. 


La energía de una señal, E, por definición es: 
co 2 
= Lio ar 


El teorema de Rayleigh, de la misma forma que el Parseval proporcionaba una 
segunda forma de calcular la potencia media de una señal periódica, proporciona otra forma 
de calcular la energía de una señal no periódica. 


E=huta= lx 


Demostración: 


ESF xPar=f xx ma = E x0|E, x pera] di 


si se tiene encuenta que una integral no es más que un sumatorio de terminos 
infenitesimales, se puede cambiar el conjugado de la integral por el conjugado del 
integrando: 


= ExojF Xx yaya ES 


cambiando de Fs las integrales: 


LE xnExorraly-Exwnmxe= Ext a 
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Este teorema es muy parecido al de Parseval, pero existe la diferencia que aquí se 
calcula la energía y allí la potencia media. Es así, porque en señales periódicas la energía es 
siempre infinita y por tanto es un dato que no aporta información, mientras, que en señales 
de energía finita, la potencia media es siempre cero y tampoco da ninguna información. 


2.5.1. — Sentido físico de |X(£)?. 


A partir del teorema de Rayleigh, se puede dar un sentido a la expresión |X(f 1% es 
la forma como está repartida la energía de una señal, en función de la frecuencia, por eso se 
suele llamar Densidad Espectral de Energía. De la misma forma cuando se hable de una 
distribución de potencia en el dominio de la frecuencia, se le llamará Densidad Espectral de 
Potencia (esta denominación se utilizará cuando se hagan cálculos sobre señales de potencia 
media finita no periódicas). 

2.6. Teoremas de la Transformada de Fourier. 


Un par transformado se suele representar por: 


x() > X(f) 


Para desarrollar y demostrar algunos de los teoremas, se seguirá siempre la 
nomenclatura: 


FlMO)=X0) y  FUXIM)=0 
2.6.1. Linealidad. 
ox(1)+ By) o aX ($)+ BY) 


Combinaciones en el dominio del tiempo se convierten en combinaciones en dominio 
de la frecuencia y viceversa. 


Demostración 


"00 ; 


[ox +prn)erra= [area pued aX) +PYO) 


2.6.2. Desplazamiento en el tiempo. 
x(1-19) 9 X(Pe 4% 


Demostración: 


UN x(t-1 Je “dt= A xAmeaAlgy= Xx (f)0 rn 
f=1t-k, 
Un retardo representa una variación lineal de la fase. 


2.6.3. Desplazamiento en frecuencia. 


NN E 
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De Í% > XLS + Lo) 
x(1)e/% > X(f fo) 


Gráficamente representa: 


x(D 


Demostración: 


fr x(1) 07404 749 qe = q x(1e dr = X(f+f,) 


«o 00 


"(aerea = |" nerd =X ff) 


o 


2.6.4. Corolario: Teorema de la Modulación. 


La propiedad anterior es muy matemática, puesto que no existen señales como las 
exponenciales complejas., pero si que existen señales tipo senos y cosenos, que con las 
fórmulas de Euler se puede descomponer en exponenciales complejas. 


Aeogt+8) , ¿2 h0001+0) 
E le YA PIO +0) 


x(t)cos(o,1 +0) = x(£) 2 


x(0 Y(D 


Y) =F(y0) =FÍx()cos(o,1)) 


Con la multiplicación del espectro de una señal por un coseno o un seno, se puede 
trasladarlo al margen de frecuencias que uno desee. * 
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2.6.5. Cambio de escala. 


Dd 
x(at) O tal x(£) 
Demostración: 


co A 00 oo jota o - joe dt 1 íí) 
=Jot ye — Ia = Pia is = 
j x(abje dt El x(atje di Sl x(t Je ps = 


t' =at 
a dt' = adt 


El valor absoluto aparece debido a que valores negativos de la constante, 
s provocarían cambios en los limites de integración, que al volver a su posición original 
cancelarían el signo de la constante. 

Una interpretación de las consecuencias de esta propiedad sobre las señales es que 
cuanto más se comprima una señal, los cambios más bruscos se harán, y esto a su vez 
provocará tener un espectro más ancho. El caso extremo sería comprimir una señal hasta 
tener una Delta de Dirac, que siendo consecuente con esta propiedad, tendría que tener un 
espectro plano, más adelante se verá que esto, evidentemente, es así. 


x(t) 


DA 


xQt) 


1/2[X(1/2)] 


Y 
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0) F(5(0) 
1 


2.6.6.  Dualidad. 


0) 9 XP) 
Xx) O x=f) 


La dualidad es una propiedad cuya utilidad se encontrará solamente en el cálculo de 
transformadas. Dicha propiedad se basa en el parecido que existe entre la Transformada de 
Fourier y su inversa. Para entenderla hay que abstraerse un poco y verla desde el punto de 
vista puramente matemático. 


Demostración: 


_ Si se calcula FHx(-f)): 
Prepera =Fxyyerray 


É f=f' 
dl = df" 


Si se compara este resultado con: 
x= nera 


desde un punto de vista puramente formal sólo hay un cambio entre las variables: 
fot 


por lo tanto, el resultado de calcular F Me f y) sería X(t). 
Si se aplica a un caso concreto 


x(t)= a) o Xx) = At, sino f4,) 


Por dualidad se tendría el par transformado: 


X()= Atysinot,) pS xf)= ak) 


En este caso, y dado que el pulso es una función par, x(f) se puede expresar como: 
x(f)= A 1) 
1 
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2.6.7. Diferenciación. 


0 o ax) 
d”x(t) "nx 
== YX) 


La propiedad de la diferenciación sólo tiene sentido utilizarla si la función derivada 
de x(t) existe. 
Demostración: 


a 7 E dd ]=f xl Lera = Pp jaen ay 


Según la definición de Transformada de Fourier, por comparación: 
dx(t 
O oia 
Derivando n veces se llega al segundo par transformado. 


2.6.8. Teorema de la convolución. 


Xy 0 O XT 
xX0y0 eo X0YY0) 


Demostración: 
La primera demostración: 


Flx(0* y (m)= FE. x(1)y(t- Dare 32% y — 


cambiando el orden de integración: 
Eolo raja xr 
ner? 


La segunda demostración: 


roo Exwomrraalrmergjea- 


cambiando el orden de integración: 
ErulEsoererala= FIX =98=Y0 YX) 


O) 


Ejemplo del teorema de la convolución. 
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x(t) 


, t 


Si se calcula directamente el espectro mediante la Transformada de Fourier: 


o : A ¿A : 
X= Pxnerta = p (+t)e“"a +/ bl 1) dr 
1 1 


. an 
- O. aprovechando que la señal es par: 


Xx) =2 (1) costara =2.; La —1,)cos(2rfiydt = 


= Af tcos(2r/1) dt -— 1l' cos(2xfdt| = 


El primer termino es una integral que hay resolver por partes, mientras que el segundo 
termino es una integral inmediata. El resultado de la integral por partes es muy conocido, 
es: 


| xcos(x)dx =cos(x) + x sen(x) 


El resultado es: 


X(f)= At, sino? (11) 


Este proceso como se ve es largo y laborioso. Sin embargo, existe una forma más 
fácil, que es utilizando el teorema de la convolución: 


As, Al) = ma) a) Sut) =v(0*v(1) 


h 


Cálculo realizado en 1.4.3 para el caso A=B y 1, = 7,. 
Como la Transformada de Fourier de un pulso ha sido calculada antes y era: 
V(S)= At, sinc(Jt) 
UN) =V IV = AN sino? (ft) 


Lo se tenía es una triángulo de amplitud A, por tanto: 
U y 
x= 2D a at sine*(4t) 
A 
2.6.9. Transformada de Fourier de S(t) y S(£). 
FÍ8(0)= Fsmentar =1 


Recuerde la propiedad de la Delta de Dirac: - 
[stoxiar=x00) 
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por tanto: 


401 


Análogamente: 
rus) =P spera =1 
1os(y) 


Evidentemente, se puede utilizar la propiedad de la linealidad para encontrar los 
pares transformados: 


AÑ) O A 
AN) 


2.6.10.  Corolario: Transformada una señal periódica. 


Partiendo del par transformado: 


AD ASÍ) 


y utilizando el teorema de desplazamiento en frecuencia: 
Aci > ASS Jo) 


Ahora, recurriendo a las fórmulas de Euler y al teorema de la linealidad se puede encontrar 
el par transformado: l 


Acos(m,1 +0) + Hera Aye "8 +/,)) 


Este par transformado es muy importante porque permite encontrar la transformada 
de Fourier de una señal periódica. Una señal periódica siempre se puede escribir en forma 
de su desarrollo en serie: 


o 
x(t)= Eee 0 =2x 


n=-0 
y teniendo en cuenta el par transformado descrito arriba: 


X(f)= Y CAS —nf0) 


n=-0 
Si, el Desarrollo en Serie de Fourier estuviera en forma de senos y coseno, se 
utilizaría el otro par transformado. 


Es importante darse cuenta que se realiza la Transformada de Fourier sobre el 
desarrollo en serie y no sobre la señal periódica inicial. 


2.6.11. Teorema integral. 


xo eo XU)JUO) 
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Demostración: 


Siendo x(t)*5(t)=x(t) y: 
9(1)= tuo) 


x(t) se puede expresar de la siguiente EA 


x(1) = «0 (uo) = OO) 


Integrando a ambos lados de la igualdad: 
xa = x(0)*u(1) 


con lo que: 


FL xoa0) = Flx(0*u(0) = Xx UY) 


Más adelante se calculará la Transformada de Fourier de un escalón, aquí de todas 
formas se utilizara la transformada del escalón porque normalmente este teorema viene 
escrito de esa forma: 

so hy 


Fr) 


con lo que el par transformado queda: 


E X(0)6 Xx 
na o 2D, Mo 


2.6.12. Ejemplos: Utilización de los teoremas de la Transformada de 
Fourier. 


Encontrar Fleu(o) 


| 2. 1 


-t -j2nA (14327 )t 
x= Lera = [rg e 


1+ ¡21f 


Encontrar Fe”) 


Pa A ¿2 oos(arfo) +27 | EE 


1+(277) +(25f)' 
Encontrar F (ent sen(£) 


Como se puede ver en la figura,la transformada que se pide calcular es de una 
función impar, y por tanto se puede utilizar la simplificación calculada para este caso. 
También se podría descomponer en dos transformadas, una para la parte negativa del eje de 
tiempos y la otra para la parte positiva, ésta última ya ha sido calculada al comienzo de este 
apartado, por tanto, solamente quedaría por calcular, la transformada para el eje negativo de 
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tiempos. Aquí se a decidido utilizar la característica de función impar para:encontrar la 


transformada propuesta. 
ie 


(—sen(2rf1) -21f | IT a 
1+ (257) | ay) 


=> 


Por el teorema del cambio de escala: 


x(at) O E x£) 


sd dE 
- jan 
JAM 


sgn(t) 


Si la función anterior se lleva al limite cuando a tiende a 0, se tiene la función signo, 


por tanto, si se lleva al limite su transformada, se tendrá la trasformada de la función sgn(t). 


A A A A 
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Encontrar ii 
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e AM 
= lim => 
Hal] 
senl(1) > = 
Encontrar FÍu(1)) 
1 1 1 1800.. 1 
3 FU+sen(o)=>3| SU+ H 1 a 
aye DD E Ls : 
a 
Encontrar la Transformada de Fourier de: 
x(t) 
A 


, t 


Esta transformada ha sido encontrada antes, ahora se calculará utilizando el teorema 
integral: 


E . |. AA = pe We =x(1) 


de esta forma se calculará primero la transformada de la derivada de x(t), que es más fácil 
de calcular, y luego, se aplicará el teorema de la integral. 


(5), Xa (1) , Xa (0500) , Aa) ] 
aan 2 Jay ] 


F(x(0) = ESO 
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EN 
xa (1)= axe lo DA at” 


1 


+2 
Xoz ()= AY ” 2 1 Xa (A) =-At, sinc(ft, Je 992 


1 


Xa (0) = At 
Xo2 (0) =-At; 
Mii a [41 sino( Jue To — Ar sino( fe 7 0P]= 
a e sinc(f/25 sen(ry,) = At,sinc? (ft,) 


que evidentemente es el mismo resultado obtenido anteriormente. 


Encontrar F* E con a >0 
(a+ ¡21f) 


Se conoce el par transformado: 
Fle"u(o) 
1 
1+ /27f 


“uno 


así mismo, se disponen del teorema del cambio de escala, que si se aplica a este par 


transformado: 
1 1 1 
Cu o ———= 
21+ 2h a+j27 


Por otra parte: 
Lane UY Xx) 
dt 


y por el teorema de la dualidad: 
d” 
df” 


Xp) o Eja) x(t) 


Si se calcula la derivada enésima de: 
1 


(a+ ¡2xf) 
d” ¡ia n(-j2m)" 


df” (a+ j2rf) (a+ j2rf)” 


Despejando la parte que interesa: 
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1 S dr | 1 | 1 
(a+ j2xf)" 27 Llar 3298) (m- Dejan)” 


de esta forma, se puede calcular la antitransformada pedida aplicando el teorema de la 
derivación y él de dualidad, sobre una transformada inversa conocida: 
1 1 po 


É (a+ 3217)" e a MT 


eur) 


Encontrar / l Ns = »D) (tren de deltas). 


n= 


Al ser periódica, encuentro el desarrollo en serie: 


00 00 1 
N ar-aT)= Y CeRmde "jp A 
n=-0 n=-00 


Calculando los coeficientes: 


Eu 


1 , 1 
EE 8lt JA Y — 
OS 


T 


Por tanto, se puede escribir el tren de deltas como: 


o) 0 1 ¡ 3 
> 1 - = > pm 2 mft 
1 —nT) pe 72 e 


n=-0w0 n=-0 '=-00 


Como se conoce el par transformado: 
A0to0tO o neió5 y — fo) 


2 E o 
ñ Sa0—00) = EAU fo) = Jo L8S 019) 
Encontrar F (l F x(1- ryeyarar) siendo x(0) e y(t), funciones pares. 


0 [E x-eynarao= [xo yaya = Lx yiar Ly ya 


Por tanto: 


a HD y) crea) 


Para calcular la constante: 
Exuyyoar =P xy = Exa = [yr = 


|f=0 
_XOYO) _ ds 
2 
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Si se sustituye la Cte. por el valor conseguido: 


A Exar-myenarar) = +. pu ) ¿ZO yy) LOYO y er 


e XUNO) 
jo 


Encontrar F(sen(ld)). 


Para encontrar esta transformada se puede descomponer en: 
- Físen(l)) =FÍsen(o) seno) = Flsento) )« Fisgn(1)) 
La transformada de la función signo ya fue encontrada, y es: 


senl) > E 


Para encontrar la transformada de un seno, se debe utilizar las fórmulas de Euler: 


EEES 


2) 


Físen(0))= ne 


Una vez encontrados los dos términos sólo queda realizar la convolución: 


dry E 2. O 2 


= 1 , 
LEIA 


2.6.13. Ejemplos propuestos. 


Encontrar la Transformada de Fourier de : 
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x(t) 


3 A 1 3 t 


Nota: Encontrar la Transfomada de Fourier de la derivada de x(t) y aplicar el teorema 
integral. También se puede descomponer la figura como resta de dos triángulos. 


Resultado: 


X(f)= 8sinc[4.f) sinc(2 1) 
Encontrar la Transformada de Fourier de: 


x(t) 
A+B 


pe 
L-2tb=t t, h+rh+2tt 
Nota: Descomponer la figura en otras más sencillas y utilizar linealidad. 
Resultado: 
X(f)=21 sinc(2 fr)e?Yo(2.4 sinc(2 ft) + B- 4) 


Encontrar la Transformada Inversa de Fourier de: 


Aja A, Aja - fa f, qq él 
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Nota: Aplicar dualidad y desplazamiento en frecuencia. 


Resultado: 


x(£) = 2asinc*(ar) - cos[27f,1) 


2.7. Función de transferencia. 


Si se recuerda la definición de respuesta impulsional, en sistemas lineales invariantes: 
la respuesta impulsional es la respuesta al sistema cuando la excitación es una Delta de 


Dirac. : 


x(1) = 6(0) > y(1) =h(0) 


y la relación entre la entrada y la salida para cualquier excitación: 


y (1) = x(0*h(t) 


Ahora, si se realiza la Transformada de Fourier queda: 


YN =XMH0) 


De esta igualdad ya se vio la expresión H(£), que es la función de transferencia del sistema o 
la respuesta en frecuencia del sistema. 


Ho)= nera 


Esta función de transferencia no es nada más que la Transformada de Fourier de la 
respuesta impulsional. Esto permitirá encontrar la respuesta de un sistema a cualquier señal 
de entrada, siempre que el sistema sea lineal e invariante y que la señal sea de potencia 
media finita o de energía finita, o sea, que tenga Transformada de Fourier. La respuesta se 
puede encontrar como la multiplicación de la respuesta en frecuencia del sistema por la 
transformada de la señal de excitación y no con la convolución, que es una operación 
complicada y laboriosa de realizar. 
A su vez, se puede presentar la relación entrada-salida de otra forma: 


Y) 
ENE 


que permite el diseño de la función de transferencia sabiendo la entrada que se dispone y la 
salida solicitada. Este cálculo puede llevar a errores, puesto que la H(£) encontrada puede 
no ser causal o fisicamente realizable. Incluso puede darse el caso que H(f) no tenga 


transformada inversa, puesto que esta función se ha hallado dividiendo otras dos y no como 
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transformada de otra función, lo que por ser la Transformada de Fourier una relación 
biunívoca aseguraría la existencia de la transformada inversa. 

De la misma forma que en señales periódicas se podía calcular la potencia media a la 
salida de una sistema L.T.I. conociendo el desarrollo en serie de la señal de entrada y 
aplicando el teorema de Parseval. Aquí se calcularán energías y se utilizará el teorema de 
Rayleigh, él cual decía: 


ES xa 


Si se aplica este teorema a la señal de salida de un sistema, y(t): 


00 2 00 2 2 
A E 
por comparación, la densidades de espectrales de energía tanto, a la entrada como a la 
salida del sistema serán: 


xl 
ro = lx ña 


No hay que olvidar que de todo este desarrollo las señales que se pueden consideran 
reales, en el sentido que existen o pueden existir en el mundo físico, son las señales 
temporales, x(t), y(t) y h(t), el dominio transformado, es decir, el dominio de la frecuencia, 
sólo es un paso intermedio, puramente matemático, que se utiliza para resolver más 
fácilmente este tipo de problemas. 


2.7.1. Formas básicas de interconectar sistemas L.T.I. 


En este apartado se calculará la función de transferencia de un sistema creado por 
conexión de varios subsistemas. Esta forma de análisis es importantes porque, en sistemas 
complicados, donde se debe realizar un estudio global, es necesario descomponer el sistema 
en subsistemas más pequeños, que puedan ser analizados o diseñados más fácilmente y 
luego interconectarlos, obteniendo el sistema completo. É 

Las formas más usuales de interconectar sistemas son: 


En cascada: 


vít) 
x(t) 0) H 16) y(t) 


Si se conectan en cascada, se tiene: 
j VO) =x(0*h(0 


y(M= o = (x(ym(0)h,(0) 


Si se realiza la transformada de Fourier: 


VA =X0MA A 


YN=/DE.O) OAMI 


como para calcular la respuesta impulsional, o la repuesta en frecuencia, la excitación tiene 
que ser una Delta de Dirac, y su transformada vale 1: 
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HR H (DA) 


Por tanto, el sistema equivalente queda: 


y HOR — 0 


Se habrá calculado por separado H,(£) y H>(£) introduciendo una 3(t) en cada 
uno de los sistema y realizando la transformada de la respuesta impulsional. 


En paralelo: 


x(t) 


VAN =XMA0) 
TS) = XD) 
YD =VIM+Z0) 


Y) = XA, +0) 


Si se sustituye X(f) por la transformada de la Delta de Dirac: 


HP = HD 


Por tanto, el sistema equivalente queda: 


xy — H(0+H(D— yo 


Realimentación: 


y). 


ZP =YWMNHO) 
A DES ARA RES RES AR 
TO =YOADN= XI +Y0MEMALN 
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A) 
NS 
DAN e 


Sustituyendo X(f) por la transformada de una Delta de Dirac: 
HI) 
HN=— ER 
1H (MAL 


De esta forma, se puede sustituir los dos bloques más el sumador, por un sistema con 
función de transferencia H(£). 

Evidentemente, todos los cálculos realizados para interconectar sistema, se pueden 
hacer en el dominio temporal, pero no son tan fáciles de resolver. 


2.7.2. Ejemplo: Cálculo de un sistema inverso. 


Se vio en el primer tema que un sistema invertible era aquel, para él que se podía 
encontrar otro sistema, tal que, conectado en cascada la respuesta impulsional total sea una 
Delta de Dirac, con un cierto retraso y con una cierta amplificación. 


h(0)*h, (6) = KS(1-1p) 
HH (S) = Ke Ho 


Por tanto, para encontrar el sistema inverso a un sistema cuya respuesta el frecuencia sea 
H(f), simplemente hay que despejar de la igualdad anterior: 
Ke Ro 


HD == Gp 


Esta operación siempre se puede realizar, lo que no significa, que el resultado 
alcanzado sea realizable fisicamente, por ejemplo, no es realizable un sistema H](f) que en 
algún punto valga infinito. 

Por otra parte no es necesario que un sistema invierta a otro en todo el margen de 
frecuencias, sino que simplemente ha de hacerlo en el margen de frecuencia en el que se va 
a utilizar es sistema, la repuesta del sistema inverso en un margen de frecuencias que no se 
utilizará es indiferente. 

El cálculo realizado en este apartado sólo es válido para sistemas lineales, puesto 
que se está utilizando conceptos como respuesta impulsional y respuesta en frequencia. Sin 
embargo, la invertibilidad de un sistema no es algo que quede restringido a sistemas lineales, 
es una propiedad mucho más general. Es posible encontrar el sistema inverso de un 
amplificador logarítmico, y este es un sistema no lineal, simplemente que el cálculo de dicho 
sistema inverso no es posible realizarlo de esta forma. 


2.73. Ejemplo: Conexión de dos sistema en cascada. 


. 


Se propone calcular la respuesta impulsional de dos sistemas como él siguiente 
conectados en cascada: 
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En el sistema de la figura se establece la condición: 
Ti>T2 
Sistema 1 


1] 
1 
J 
1 
1 
1 
1 
1 
| 
¡x(t) 
1 
1 
1 
] 
1 
I 
1 


Primeramente se debe analizar un sistema por separado y después se interconectarán 
los dos sistemas. ¡ 
Se calcula la respuesta impulsional, excitando el sistema 1 con una Delta de Dirac. 


v(t) = (1-1,)- At TZ) 


2 TAE) 
y) = 0 =Í[5(—7)-5(0-7 Jai =u(1-7)-u01=2)=- + 


Si ahora, se encuentra la respuesta impulsional: 
T+T, 


HS)= Fo) =A7-2)sind y(1-z) Je? 


Por último, se conectan los dos sistemas en cascada: 


x(t Sistema 1 z(t 
(0 Sistema y (y 


20) = (OO), z, MOMO 


en términos de frecuencia: 
0% -j: + 
EN = HAN = (1-7) sino [7-7 je rra 
y se realiza la antitransformada de Fourier: 


nos AN—=2 
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2.7.4. Ejemplo: Cálculo de la respuesta en frecuencia de un amplificador. 


Se dispone de un amplificador de audio, del cual se quiere saber su respuesta en 
frecuencia. Para lo que se realiza el siguiente montaje. 


Analizador 
Generador Amplificad. de 
espectros 


Dada la imposibilidad física de introducirle una Delta de Dirac, para que después el 
Analizador de Espectros realice su Transformada de Fourier, se decide utilizar un generador 


que proporciona la siguiente señal: 
É 
x(t) = al) 
T 


El generador permite variar la anchura del pulso, es decir, 7, y la amplitud del pulso A. 
Evidentemente, utilizar un pulso y no una Delta de Dirac producirá un error en la 

respuesta en frecuencia medida. ¿Cuánto tiene de valer T, para que el error cometido a 

cualquiera de las frecuencias donde se realiza la medida sea inferior a un 1%, si además, se 

quiere que la respuesta del sistema a la continua sea exactamente la medida? 

Nota: El Analizador de Espectros es un instrumento que calcula y representa la 

Transformada de Fourier de una señal. 


Solución: 


Se supone que el amplificador tiene una respuesta en frecuencia H(£), por tanto, la 
señal que medirá el Analizador de Espectros será: 
YN =HNXA) 
Y(f) = E) Arsind(fr) 


El resultado obtenido se está multiplicando la respuesta en frecuencia del 
amplificador por una función sinc, que no pondera por igual todas las frecuencias. Si la 
respuesta en frecuencia del amplificador es constante a todas las frecuencias, cosa deseable, 
se obtendrá la multiplicación de las dos funciones superpuestas en la siguiente gráfica. 


De un primer análisis, se puede observar que el valor de 1 tiene que ser lo 
suficientemente pequeño como para que el primer cero de la función sinc no este en el 


Página 45 


¡Análisis de Fourier : Apuntes procesado de señal 
-—==z==---5--AA=2=2=1 683 yáxy” 0 O y putes procesado de señal 


margen de frecuencias donde se realizará la medida, puesto que si está dentro: el error 
cometido será demasiado grande, sería del 100%. 
Por otra parte, el error será mayor cuanto más alejado del origen de coordenadas se 
esté. Por consiguiente, la condición de error máximo se impondrá en la máxima frecuencia a 
medir, como el amplificador es de audio, la frecuencia máxima son 20K Hz. 
Si se define el error como un error relativo: 
medida correcta - medida con error 


medida correcta 


error relativo = 


Conocida la medida correcta H(£) y la real Y(£), solamente hay que imponer la 
condición que el error sea menor de un 1%: 
H(0KHz) - Atsinc(20XH71)H(20KHz) 
H(20KHz) 


100=1 


Imponiendo la otra condición, que el valor de continua sea correcto: 
H(0Hz)- A WH(0KHz) _ 


0>47T=1 
H(OKHz) a 
queda: 
sen(2000077) 0:00 
2000077 


Si se resuelve esta ecuación, por cálculo numérico, o recurriendo a la tabla de 
valores de la función sinc, se llega al resultado: 
71 3.12 10 seg. 
A=320KV. 


Ahora, se supone que el generador no es capaz de dar esos valores, ¿cómo se 
modificaría el montaje, si no se quiere cometer ningún error? . 


Solución: 


Si no se modifica el montaje se está midiendo: 


YAH=HAMX0) 


mientras que lo que se quiere medir es: 


YAN=H 


Por tanto, se tendrá que introducir un sistema que permita compensar los errores 
que introduce la transformada del pulso. La solución es introducir un sistema de tenga una 
respuesta en frecuencia inversa a la del pulso, para tener: 


Yy)= ENIDGF =H(f) 


Si se tiene en cuenta la transformada del pulso se tiene que el sistema tiene que ser: 
el de a 


X(f) 7 Arsincl ft) 
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Este filtro no es realizable, puesto que tiene que valer infinito en valores de 
frecuencia iguales a los múltiplos del inverso de 7. Pero como la medida tiene que realizarse 
para valores de frecuencia menores de 20KHz y se ha elegido 1 para que el primer cero esté 
fuere de la banda de 0-20KHz, se puede realizar un sistema, con la respuesta en frecuencia 
igual al inverso de la función sinc, pero solamente hasta los 20K Hz. Asi, el filtro que hay de 
colocar antes del Analizador de Espectros es S: ejemplo: 


el ED O Y. a 
us — Arsinc At) Her) 


X(F)0KHz 
Véase que se aplica la condición de sistema inverso de forma, podría decirse 
práctica, es decir, solamente en el margen de frecuencias que interesa, tal como se comentó 
en el apartado 2.7.2, 


2.8. Densidad espectral y correlación. ' 


2.8.1. Producto escalar de dos funciones. 


El producto escalar de dos funciones xx(t) e y(t), del mismo tipo, energía o potencia 
media finita, se define como: 


Exa) y" (1)dt para señales de energía finita 


DI e 
lim FL x(1) y" ((Jdt para señales de potencia media finita 


En señales periódicas, se simplifica siendo: 
1 s 
<> x0y 04 
De la definición de producto escalar se puede sacar un consecuencia directa: el 


producto escalar de una función consigo mismo es igual a la energía o potencia media de la 
señal. j 


<x(0,x(0) >= F x(0x" (Ndt = Pol ar señal deenergía finita 
1 o ; 1 o . 
<x(0, (0) >= limA) x(0x"(Ndt = lim = x(t) * di señal de otencia media finita 
T>0 T “o T>w0 T “o Pp 


Pasando a ver la interpretación que se puede dar al producto escalar, para ello 
conviene apoyarse en la desigualdad de Schwartz, la cual dice: 


E x(1) y mal s Pleojaf| yojar 


Para que se cumpla el signo igual la señales x(t) e e deben cumplir: 
y (1) = ax(t) 


Por tanto, el producto escalar de dos señales está acotado por: 
0<< x(0), y(1) >< Foxipa 
En cierta forma el producto escalar está midiendo el parecido entre dos señales. 
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Si se tiene presente el producto escalar de dos vectores, se puede dar otra 
interpretación del producto escalar de dos funciones, el producto escalar da un valor 
proporcional a cuántas veces está contenida una función en la otra. Dicho de otra forma, 
cuál es la "proyección de una función sobre la otra", 

Si el producto escalar es igual a cero se dice que las señales son ortogonales. 
Analizando los motivos por los cuales dos señales pueden ser ortogonales, se ve que hay 
dos motivos evidentes: : ; A 

- Dos señales son ortogonales sí tiene simetría opuesta, es decir, una par y la otra 
impar, 

- Dos señales son ortogonales sí no se solapan nunca en el tiempo o en frecuencia. 

Hay señales que no cumplen ninguna de estás condiciones y son ortogonales, por 
ejemplo dos:cosenos de diferente frecuencia, pero en los dos casos mencionados siempre lo 
son y son fáciles de deducir a partir de la representación gráfica. 


2.8.2. Correlación cruzada. 


La correlación cruzada es un generalización del producto escalar de dos funciones. 
No mide sólo la similitud entre dos señales, sino que además, da información sobre posible 
similitudes con un cierto desplazamiento temporal entre ellas y mide este desplazamiento. 

Haciendo un símil geométrico daría información sobre el parecido de figuras 
situadas en el espacio con un patrón y daría su posición. 

En general la correlación de señales de energía finita está definida como: 


Ro (1) = y x(1) y" (1 vdt =< x(0), y(t 1) >= x(1)* y" (-1) 
RO =P yx aa (0,x(1=1) >= peje (er) 


Como se puede observar, la correlación cruzada no tiene la propiedad conmutativa, es 
decir, en general: 
Ry (0 +Rya (7) 


Algunos autores definen la correlación cruzada de dos señales de la siguiente forma: 


Rm) = Fco +1)y" (Nat 


El significado y los valores son los mismos en los dos casos, es lo mismo retrasar y(t) 
respecto a x(t), que adelantar x(t) respecto a y(t). Ñ 
Cuando las dos señales son la misma la correlación cruzada recibe un nombre 
especial, se denomina autocorrelación: 
RO) =R, (0) = Pxore -vdt 
Hasta aquí se ha definido la correlación en señales de energía finita. En señales de 
potencia media finita, la correlación de define de forma muy parecida. 


Ya (0) = tm Ex (—ndt=<x(0,yU=1)> 


ya = pm yO (A =< 0, x1=1)> 


. 


Para el caso en que las señales sean periódicas la expresión se simplifica, siendo: 


_EFO AA <á<—<>+m 
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“y lt) = +, x(1) y" (tv) dt =< x(0), (11) > 


Todas las observaciones realizadas sobre la correlación de dos señales de energía 
finita son perfectamente aplicables a las señales de potencia media finita. 

Hay una característica especialmente importante en la autocorrelación de señales, 
tanto de energía, como de potencia media finita, si bien, éstas últimas sólo para las que sean 
periódicas. Se vio que la densidad de energía/potencia media (señales periódicas) de una 
señal venía dada por la expresión IX(f y. Una vez conocida la autocorrelación, esta misma 
expresión se puede calcular como la Transformada de Fourier de la autocorrelación. Se verá 
.. Primero para señales de energía finita y luego para señales de potencia media finita 

periódicas. ñ 
Si se llama a la densidad espectral de energía como Sxx(f), se tiene: 


RiO=R (0 = xr nar E x(0) Fx peta] a és 
=' xo] $ XP mera lar = Fx x(De aora = 


=P xy) xy yema Plena =Fs.era = FS.) 


: Por tanto, queda demostrado el par transformado: 
Ra) O Sa (S) 


En el caso de señales de potencia media finita periódicas: 


Y, (0) = y, (1) = 5 f a x(1)x" (1- v)dt= +, , x()] Ec a = 


= +E 2) Ejerce t= Ed: E pera ams = 
o po 


CS, 


0 00 2 
A id JNOT _ JnoyT 
- Pace loe 


n=-00 n=-00 
Si se realiza la Transformada de Fourier a ambos extremos del desarrollo: 


Gat1)= 20 Pas =) 


Por otra parte, se tiene que la transformada de una señal periódica era: 


x= Eca1=05) =[Xx0 Y = Eca(s=05,) = 
= Ecals 01) Ecrls 04) = le. af) 


2 


o 


C, 


De esta forma queda demostrado para señales de potencia media finita periódicas. 
En el caso que las señales no sean periódicas puede que no coincidan la transformada de la 
autocorrelación y el módulo al cuadrado de la Transformada de Fourier de x(t). En el caso 
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que no coincidan, la densidad espectral de potencia media será la transformada de la 
autocorrelación y no el módulo al cuadrado de la transformada de la señal. 


2.8.3. Propiedades de la correlación. 


DRY RED (Si (Y) 


Sx(DY(0 Ru) =Ru(-0>S.(f) es real 

Sí, además, x(t) es real, Rxx(t) es real y par. 

Todas estas afirmaciones se pueden demostrar fácilmente sustituyendo las 
correlaciones por su expresión integral. 


2) Ryx(t) no depende del origen de tiempos solamente del retardo de una señal respecto a 
otra. 


3) Rxx(t) es máxima en el origen O] < R.(0) 
Por otra parte, la autocorrelación en el origen es la energía de la señal. 


Es. = E s.Mera|_ 2.0) 


4) En sistema L.T.I, las correlaciones de las señales de la entrada y la salida están 


relacionadas: 


y0=x0:0=Í xy. 


Las distintas correlaciones de las señales de entrada y salida son: 
Ro (1) =x(1)x" (1) 
Ry kt) = yA (E) = (A (a) = Ra) 
Ry) =x(0) y (=1) = (0) (ae (a) = Ro (eh (—x) 
RO = AY a) = AAA (e (a) = Re) A (a) (a) 


Con lo que en el dominio de la frecuencia queda: 


Se (1)= 80 (HL) 
SyU)=S NE) 
SUS NIE? 


Para concluir, es importante tener en cuenta que a diferencia del espectro. de una 
señal, su densidad de energía, o en su caso, la de potencia media finita, no es única para 
cada señal, puede haber dos señales diferentes que tengan la misma función de densidad. 

La densidad espectral no da tanta información como el espectro de una señal, pero a 
veces no es posible calcular un espectro y es lo único que se puede calcular, por ejemplo, 
las señales aleatorias. . 


á ÁA A — _—=—=+4+=+141444 
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